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ABSTRACT 



Aims. The recent introduction of high-resolution/large spectral-range spectrographs has provided the opportunity to investigate the 
chemical composition of classical Cepheids in detail. This paper focusses on new abundance determinations for iron and 6 light metals 
(O, Na, Mg, Al, Si, Ca) in 30 Galactic Cepheids. We also give a new estimate of the Galactic radial abundance gradient. 
Methods. The stellar effective temperatures were determined using the method of line depth ratios, and the surface gravity and the 
microturbulent velocity v, by imposing the ionization balance between Fe I and Fe II with the help of curves of growth. Abundances 
were calculated with classical LTE atmosphere models. 

Results. Abundances were obtained with RMS accuracies of about 0.05-0.10 dex for Fe, and 0.05-0.20 dex for the other elements. 
Cepheids in our sample have solar-like abundances, and current measurements agree quite well with previous determinations. We 
computed "single zone" Galactic radial abundance gradients for the 8-12 kpc region and found a slope for iron of -0.061 dex kpc" 1 . 

Key words. Stars: abundances - Stars: supergiants - Galaxy: abundances - Galaxy: evolution 



1. Introduction 

Galactic abundance gradients are one of the most important 
inputs for the numerical models of the chemical evolution of 
Galactic disks. Different stellar tracers have been used to es- 
timate the gradient, namely planetary nebulae, giant stars, old 
open clusters, HII regions, and young, Z?-type main sequence 
stars. If the existence of the gradient is now widely accepted, its 
precise value still needs to be established. 

As a matter of fact, empirical estimates give quite steep 
gradients for ir on ranging f rom - 0.05 dex kpc 1 to -0.07 dex 
kpc" 1 : see. e.g.jFriel et al. I d2002l) who found a slope of -0.06 
al. I (12005 



dex kpc 1 ; iKovtvukh et al. (2005), who obtained a gradient of 



2006), who estimated a slope 



-0.06 dex kpc 1 ; and lLuck et al 

of -0.068 dex kpc -1 using classical Cepheids as chemical trac 
ers. However, much s hallower slope s have also been estimated 
-0.023 dex kpc' 1 by iTwarog et al.l d 1997b . -0.044 dex kpc 



by lAndrievskv et al. I d2004l) . And in a few cases the gradien t 
was not detected: -0.01 7 dex kpc" 1 by iNeese & Yosstdl988l) 
and -0.003 dex kpc" 1 bv lKilian-Montenbruck et al. I d 19941) . The 
broad range of empirical radial gradients also applies to the other 
light metals adopted in current investigations. Thorough reviews 
of gradient values for many ele ments can be found for example 
in lAndrievskv et al. I d2002al) and lChiappini et all d2001h . 

Beyond an expected scatter due to measurement errors, these 
discrepancies between steep and shallow slopes of the Galactic 



Send offprint requests to: B. Lemasle 

* Based on observations made with the 1.52m ESO Telescope at La 
Silla, Chile. 



radial abundance gradient might trace its evolution as a function 
of stellar ages. This working hypothesis has been su ggested in 
several studies dFriel et al. 1120021: iMaciel et ai7ll2003l) based ei- 
ther on planetary nebulae or on open clusters. From this point of 
view, Cepheids are young Galactic objects that trace the present- 
day abundance gradient, and only a few eleme nts (C, N, Na, 
and perhaps Mg and Al); dKovtvukh et al. 1 11996) may have their 
abundances modified during the life of the Cepheids, while all 
the others should trace the abundance pattern as it was at the 
birth of the star. 

Compared to other methods, the use of classical Cepheids 
offerss many advantages because these stars 

- are excellent distance indicators; 

- are luminous enough and provide the opportunity of deter- 
mining the gradient over a significant fraction of the Galactic 
disk; 

- allow the measurements of many elements that cannot be in- 
vestigated with other stellar tracers; 

- show spectra with many absorption lines and, in turn, reliable 
abundance values. 

An extension of the available data is mandatory for con- 
straining the slope of the Galactic abundance gradients on a more 
quantitative basis. This paper is the first in a series focussed on 
accurate and homogeneous abundance measurements for a good 
sample of Galactic Cepheids. The sample includes objects that 
are in common with previous studies, but current spectra typi- 
cally present a higher spectral resolution and/or a wider spectral 
range. 
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2. Observations 

The observations were performed in 2002 at the 1,52m ESO 
telescope in La Silla. We collected 38 high-resolution spectra 
of Galactic Cepheids with the fibre-fed Extended Range Optical 
Spectrograph (FEROS), with a resolving power of 48,000 and a 
spectral range spanning from 370 to 920 nm. The spectra were 
reduced using the FEROS package within MIDAS. The S/N ra- 
tio at 600 nm ranges from 40 to 136 (see TableQ}. 

3. Methodology 

The determination of an accurate temperature is a critical point 
in the abundance determination. In this investigation, the temper- 
ature was derived spectroscopically and does not rely on photo- 
metric indices. 



2. The quality of each fit is estimated by calculating the^- 2 . 

3. A new "generation" of Gaussians is calculated from the 20 
best fits after adding random modifications to the initial set 
of parameters values ("mutations"). 

4. The new set replaces the old one, and its accuracy is esti- 
mated by once again using a x 2 evaluation. 

5. The process is iterated 100 to 200 times until the conver- 
gence to the best Gaussian fit is attained. 

In the current investigation we limit our analysis to weak, non 
blended lines, all of which could be fitted by a Gaussian. The 
strong lines were neglected. Some weak but asymmetric lines 
were also neglected. This selection reduced the number of lines 
really used to 30-90 for Fel, 5-10 for Fell, 1-3 for O, Na, Mg 
and Al, 1-10 for Ca, and 10-20 for Si. 



3.1. Line list 



For ir on, we used the line list proposed by iMottini et al. I 
(2007), which contains 246 Fe I lines and 17 Fe II lines 
covering a wide spectral range corresponding to the FEROS 
one. These lines are a compilation of the l ine lists from 
Cle mentini et"aTI (1 19951). iFrv & Carney I (1 19971) . iKiss & Vinkol 
(2000) and lAndrievskv et al. I d2002al) . Thev were completed 
with a sel ection of lines from VALD (Vienne Ato mic Line 
Database) dKupka et al. ll999l:lRvabchikova et al. 1 1999b . picked 
for effective temperatures corresponding to Cepheids: 4500- 
6500 K. 

Mottini et al. devised the following approach to removing 
blended lines. 

- All the lines stored in VALD between 4800 and 7900 A and 
corresponding to Cepheids typical parameters have been 
compiled. 

- This compilation was over-plotted onto the iron lines, with 
special attention given to the lines located within + 3 A of an 
iron line. 

- If the these lines contributed to more than 5 % of the iron 
line strength, then the iron line was neglected. 

This test was performed on three observed spectra for Cepheids 
with effective temperatures about 4500 K, 5500 K, and 6500 K. 
It might remain a few slightly blended lines, since a few weak 
lines may not be listed in VALD, but their overall contribution is 
negligible. 

The same method was adopted to clean the line list of alpha 
elements, which are much less numerous: 4 lines for All, 9 lines 
for Mgl, 4 lines for Nal, 2 lines for 01, 20 lines for Cal and 24 
lines for S il. For all the lines, the physical properties (oscillator 
strength, excitation potential) listed in VALD were adopted. 



3.2. Equivalent widths 

The measurement of equivalent widths (EW) was performed 
with a semi-interactive code (fitline) developed by R Francois. 
This c ode is based on genetic algorithms from Charbonneau 
d 1995b . and it mimics the way genetic mutations affect DNA, 
driving the evolution of species. Lines are fitted by a Gaussian 
defined by four parameters: central wavelength, width and depth 
of the line, and the continuum value. Line after line, the algo- 
rithm runs in the following way: 

1 . An initial set of Gaussians is computed, giving random val- 
ues to the four parameters. 



3.3. Temperature determination 

We obtained the stellar effective tempera ture T e ff by using the 
metho d of line depth ratios, described in iKovtyukh & Gorloval 
(2000). This method only relies on spectroscopy and, therefore, 
is independent of interstellar reddening. The key advantage to 
this approach is that it uses weak, neutral metallic lines appro- 
priately selected: 

- The excitation potentials of line depth ratios must differ as 
much as possible for a good estimation of T e /f, but lines 
should be relatively weak to limit the dependence on micro- 
turbulent velocity. 

- Pairs of lines must be close enough to eliminate errors in 
continuum placement. 

- Ionized lines are excluded because they are too sensitive to 
the surface gravity. 

- The line depth ratios are based on weak lines, and thus metal- 
licity effects on temperature determination are strongly mod- 
erated. It is worth noting that two lines located on different 
parts of the curve of growth do not have the same behaviour 
when metallicity changes. 



IKovtyukh & Gorlova I d2000h proposed 32 analytical relations 
for determining T e ff from the line depth ratios. After measuring 
the line depths and calculating subsequent ratios, we obtained 32 
temperatures whose mean value gave us the stellar temperatures 
listed in Table|2] According to the quoted authors, the uncertain- 
ties on absolute temperatures are of the order of 80 K, and un- 
certainties on relative temperatures inside a sample are smaller 
and of the order of 10-15 K. 



3.4. Atmospheric parameters 

The surface gravity, log g, and the microturbulent velocity v f are 
determined by imposing the ionization balance between Fe I and 
Fe II with the help of the curves of growth. Lines of the same 
element in different ionization states should give the same abun- 
dance value. Iterations on the values of logg and v, are repeated 
until Fe I and Fe II adjust to the same curve of growth. Figures[T] 
and|2]show the curves of growth for Fe I and Fe II in AP Pup. An 
independent support to the T e jf value is given by the evidence 
that both high and low^„ lines properly fit the curve of growth. 
Once abundances were calculated, atmospheric parameters were 
checked and in particular, the Fe I abundances depend neither 
on line strength nor on excitation potential. The atmospheric pa- 
rameters for the target Cepheids are listed in Table [2] 
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Table 1. Observations: date, exposure time, and S/N at 600 nm. 



Object 


date 


exposure time (s) 


S/N 


AD Pup 


02/01/01 


600 


65 


AH Vel 


29/12/00 


300 


1 1 1 




31/12/00 


100 


88 


APPup 


30/12/00 


500 


119 




31/12/00 


300 


105 




01/01/01 


300 


114 




02/01/01 


120 


53 


AQ Car 


02/01/01 


240 


76 


AQPup 


01/01/01 


500 


84 


AT Pup 


30/12/00 


600 


73 


AX Vel 


31/12/00 


400 


94 


BG Vel 


30/12/00 


600 


100 


BNPup 


02/01/01 


600 


62 


DR Vel 


01/01/01 


600 


95 


lCar 


29/12/00 


180 


42 




30/12/00 


120 


41 




31/12/00 


80 


40 




01/01/01 


120 


37 


MYPup 


31/12/00 


100 


94 


RS Pup 


31/12/00 


240 


48 


RYCMa 


29/12/00 


240 


84 




30/12/00 


600 


67 


RY Vel 


29/12/00 


600 


69 


RZCMa 


02/01/01 


600 


77 


RZ Vel 


01/01/01 


300 


42 


ST Vel 


02/01/01 


600 


80 


SW Vel 


31/12/00 


500 


153 


SX Vel 


31/12/00 


400 


82 


T Vel 


31/12/00 


400 


106 


TWCMa 


02/01/01 


600 


66 


UXCar 


31/12/00 


400 


105 


V397 Car 


02/01/01 


300 


86 


VCar 


29/12/00 


600 


78 


V Vel 


31/12/00 


300 


136 


VXPup 


02/01/00 


300 


77 


VYCar 


31/12/00 


240 


130 


VZ Pup 


01/01/01 


600 


55 


WXPup 


01/01/01 


300 


53 



3.5. Abundance determination 

Abundances are t hen ca l culate d with the help of the atmo- 
sp heric model of Spite (1 19671) based on the grid of models 



of iGustafsson et al. I d 19751) . These models are based on the 
following assumptions: plane-parallel stratification, hydrostatic 
equilibrium, and LTE. The abundance determination codes ad- 
just abundances until matching the observed equivalent width. 
We adopted the scaled -solar chemical abundances suggested by 
iGrevesse et al. I (1 19961) . 

For every star in the sample, the final abundance was ob- 
tained by calculating the mean value of the abundances found 
for each line. The use of a substantial number of "clean" lines 
drives to very small error bars, which are less than 0.10 dex for 
most of them. Oxygen and aluminum abundances were checked 
with synthetic spectra, to account for the occurrence of blended 
lines. In most cases, the influence of the blended lines was negli- 
gible. In a few cases, abundances were modified, but the change 
never exceeded 0. 1 dex. 



3.6. Uncertainties in abundance determinations 

The first source of uncertainty comes from the extraction of data 
from the spectra: the equivalent width and the continuum place- 
ment (which influences the measure of the equivalent width) 
must be precisely determined. To minimize the errors, we limited 
our analysis to symmetric and unblended lines and neglected the 
most crowded parts of the spectra in selecting the line list. We 
only adopted weak or relatively weak lines (EW < 120 mA) and 
assumed they could be fitted by Gaussians. 



The second crucial point is to have precise temperature deter- 
minations, because this intrinsic parameter strongly affects the 
line strength . We used, as mention e d abov e, the calibration re- 
lations from iKovtv ukh & G orlova (2000). According to these 
authors, the uncertainties on absolute temperatures are at most 
~ 80 K. Since current spectra have been secured with a spectral 
resolution of only 48000, the intrinsic dispersion in temperature 
values is on average « 80 - 130 K. 

In order to constrain how the uncertainties affect abundance 
determinations, we finally computed abundances by adopting 
over or under-estimated values of the atmospheric parameters. 
The test was carried out with AP Pup, one of our spectra with 
the lowest S/N (S/N=53). As expected, only errors on tempera- 
ture had a noticeable effect: an error of 100 K on the effective 
temperature implies a difference of about 0. 1 dex on abundances 
values, while errors on surface gravity of Alogg = +0.3 and on 
turbultent velocity of Av, + 0.5 km/sec cause much smaller dif- 
ferences, namely 0.02 and 0.05 dex. 

Finally, the LTE assumption in the atmosphere models might 
not be appropriate for supergiant stars, hence the adopted stellar 
atmosphere models could provide abundance determinations af- 
fected by deceptive systematic errors. This effect should be min- 
imized, however, by the fact t hat we only adopted w eak or rel- 
atively weak lines. Moreover. iFrv & Carney I Q997) show that 
a canonical spectroscopic approach, using classical LTE atmo- 
sphere models, gives reliable abundances for Cepheids. In par- 
ticular, they noted that LTE analysis of dwarfs and Cepheids lo- 
cated in the same open c luster s give, within the errors, the same 
metallicitv. lYong et al. I C2006) used the same method (classical 
LTE models) and examined the difference [Til/Fe] - [Till/Fe] to 
test the influence of NLTE effects. The discrepancy between the 
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two different approaches is -0.07 dex and lies within the mea- 
surement uncertainties, thus suggesting that the surface gravities 
derived from ionization equilibrium of Fe spectra are robust. 

3.7. Distances 

Distance c alculations are based o n NIR (J, H, K-bwA) photom- 
etry from lLanev & Stobie I (119941) . For a few objects we adopted 
NIR photometry from the 2MASS catalogue. The mean mag- 
nitude of these obj ects was estimated using the template light 
curves provided by Soszvris ki et al. I {2005), together with the 
V-band amplitude and the epoch of maximum available in the 
literature. Laney & Stobie (1994) magnitudes have already been 
un-reddened, while for un-reddening the new data, we first com- 
puted the relative ab sorption in the V-ban d, Ay = 3.1 E(B-V), 
and then according to lCardelli et al. I dl989h the absorption in the 
J, H, and ^-band, namely Ay = 0.28 Ay, Ah = 0.19 Ay, and A* 
= 0.1 14 Ay. For individual Cepheids we adopted the red dening 
values given in the Fernie's database dFernie et al. 1119951) . 

The Cepheid absolute magnitudes were computed using 
the NIR perio d -lum inosity (PL) relations recently provided by 
iPersson et al. (2004), together with an LMC distance modu- 
lus of 18.50. Pulsa tion periods are from the Fernie's database 
dFernie et al. 1 1 19951) . The period of Cepheids pulsating in the 
first overtone were fundamentalized, while we adopted the fun- 
damental period for the double-mode pulsators pulsating si- 
multaneously in the fund amental and in the fir st overtone. The 
PL relations provided by IPersson et al. I d2004l) are based on a 
large and homogeneous sample of LMC Cepheids, but it can 
be applied to Galactic Cepheids. Indeed, at NIR wavelengths, 
the slopes and the zero -points of the PL relation depend very 
weakly on metallic ity dBono et all 119991: IPersson et al7ll20o4 
iMottini et ani2007[) . Moreover, the NIR PL relations present an 
intrinsic dispersion that is significantly smaller than the disper- 
sion of the optical PL relations. The difference is due to the fact 
that the former ones are, at fixed period, marginally affected by 
the intrinsic width in temperature of the instability strip. Cepheid 
distances based on NIR PL relations are more accurate, when 
compared with optical ones, because they are also marginally 
affected by the uncertainties in reddening estimates. 

The heliocentric distance d can now be computed: 

d = io - 2 K-^+ 5 ) , where x = J, H, K. 

The heliocentric distance adopted in the following is the mean 
value between the three distances in J, H, and K— bands to limit 
random errors. The Galactocentric distance is finally calculated 
with the classical formula: 



Rr 



[ R C hel + ( d cosb ) 2 " 2 R C,hel d COS b COS 1 ] 1/2 



where Rcj te i is the Galactocentric distance of the Sun, d the he- 
liocentric distance of the Cepheid, while / and b ar e the Galactic 
longi tude and latitude, listed in the Fernie database dFernie et al. I 
1995 ). We assumed a galactocentric distance of 8.5 kpc f or the 
sun dKerr & Lvnden-Belllll986tlFeast & Whitelocklfl997t) . 

The results relative to current Cepheid sample are given in 
Table [3] Distances based on NIR photometry are compared with 
distances based on th e catalogue of light-cur ve parameters and 
distances provided by Berdni kov et al. I d2QQ0t) . In this catalogue, 
the heliocentric distance was calculated using a period-colour 
relation and photometric measurements in the optical bands. 
Results show good agreement for Cepheids near the solar cir- 
cle and are slightly more discrepant toward the outer disk, where 
the differences can reach 0.3 kpc and even 0.6 kpc (5%) for AD 
Pup, the most distant Cepheid in our sample. 



4. Cepheid metal abundances 

Exhaustive studies of Cepheid chemical compositions are re- 
cent after rapidly developing since the introduction of high- 
resolution, wide spectral-range spectrographs. Precise values of 
Cepheid abundances are necessary for constraining the metallic- 
ity dependence of optical and NIR PL relations and for deter- 
mining Galactic abundance gradients. This investigation deals 
with accurate abundances of iron and several light metals for 
a good sample of Galactic Cepheids. The results are listed in 
Tables 4 and 5, and for each target we give the abundances and 
the standard deviations (when more than 1 line was adopted). 
These tables also include previous abundance measurements 
when available. The literature abundances come mainly from 
the systematic investigation by Andrievsky and collaborators 
who published accurate chemical abundance s for a large sam- 
ple of G alac tic Cepheids in a series of papers: [An drievskv et al. I 
(l2002allbllcT) . iLuck et al. I d2003l) . lAndrievskv et al. I d2004l) . 

Current measurements, when compared with the quoted 
abundance analyses, present on average slightly better standard 
deviations. We typically used a smaller number of lines, but their 
fits have been interactively checked one-by-one. The standard 
deviations are on average smaller than 0.10 dex for Fe and Na, 
and they range from 0.05 to 0.20 dex for O, Al, Mg, S i and Ca. 
As a whole, we found that Cepheids in the solar neighbourhod 
have solar-like abundances both for iron and a-elements, 
as expected fro m previous abundance analyses on Cepheids 
and supergiants dLuck & Bond I [19891: ILuck & Lambert 11 1 992; 
[ Luck & Wepferll 19951: lAndrievskv et al. Il2002allblld [Luck et al. I 
l2Q03h lAndrievskv et al. 112004 iKovtvukh et al. II2005I) . 



4.1. Iron 

For iron (see TableHJi, our results show excellent agreement with 
Andrievsky's, and the discrepancies are typically smaller than 
0. 10 dex. Larger discrepancies 0.2 - 0.3 dex) have only been 
found for two targets, namely TW CMa and VZ Pup. Similar dis- 
crepancies have also been found for other elements. These dif- 
ferences might be due either to different signal-to-noise spectra 
or to the adopted atmospheric parameters. However, the curves 
of growth for these two targets and related tests appear to be very 
consistent. Figures [3] |4] and [5] |6] show the results for TW CMa 
and VZ Pup, respectively. A glance at the data plotted in these 
figures shows good agreement between the predicted curves of 
growth and current measurements for both Fel and Fell lines. 
These distant stars deserve further investigations, as they are 
quite distant from the Galactic centre, and therefore very use- 
ful to firmly constrain the metallicity gradient in the outer zones 
of the Galaxy. 



4.2. Oxygen 

The oxyg en is slightly deficient, as expected after the investiga - 
tions bv lLuck& Lambert I dl98lUl985l). lLuck & Bond I dl989h . 



and the theoretical work of ISchaller et al. I d 19921) . Current re- 
sults show very good agr eement with O abundances provided by 
lAndrievskv et al. I d2004), but with the exception of SX Vel and 
once again VZ Pup . The ratio [O/Fel is ve ry c lose to the abun- 
dance s provided bv lLuck & Wepfer] d 1995b and lKovtvukh et al.1 
(120051) . However, oxygen abundances show a higher dispersion 
when compared with other elements because the lines we used 
(A = 6300 A and 6363 A) are weak and blended. In order to re- 
duce ther intrinsic dispersion we also adopted synthetic spectra, 
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Object 


P 


phase 


T.// 


logg 


V, 


[Fe/H] 




d 




K 




km/s 


dex 


AD Pup 


13.594 


0.406 


5290±146 


1.0 


2.4 


-0.4 


AH Vel 


4.227171 


0.408 


6080± 80 


1.5 


3.5 


-0.1 




4.227171 


0.887 


6483±168 


1.6 


3.2 


-0.1 


APPup 


5.084274 


0.733 


6180± 84 


2.2 


4.0 


-0.2 




5.084274 


0.926 


6148± 93 


1.7 


3.5 


-0.2 




5.084274 


0.130 


5723± 87 


1.6 


3.7 


-0.2 




5.084274 


0.329 


5557± 86 


1.7 


3.5 


-0.2 


AQ Car 


9.76896 


0.528 


5374±104 


0.5 


3.5 


-0.1 


AQPup 


30.104 


0.309 


5304±133 


0.3 


3.3 


-0.4 


AT Pup 


6.664879 


0.594 


5348±144 


0.7 


3.2 


-0.4 


AX Vel 


3.6731 


0.598 


5995±161 


1.8 


3.2 


-0.3 


BG Vel 


6.923655 


0.895 


5957±111 


1.5 


3.5 


-0.2 


BNPup 


13.6731 


0.116 


5615±81 


1.0 


3.0 


-0.1 


DR Vel 


11.1993 


0.179 


5183± 87 


0.5 


3.0 


-0.2 


ICar 


35.551341 


0.234 


4885± 87 


0.7 


3.5 


0.1 




35.551341 


0.262 


4848± 91 


0.5 


3.5 


0.1 




35.551341 


0.290 


4900±105 


0.5 


3.5 


0.1 




35.551341 


0.318 


4821± 87 


0.6 


3.2 


0.1 


MY Pup 


5.695309 


0.904 


6114±110 


1.4 


2.8 


-0.25 


RS Pup 


41.3876 


0.826 


5056±110 


0.2 


3.2 


-0.1 


RYCMa 


4.67825 


0.257 


5666±131 


1.2 


3.0 


-0.2 




4.67825 


0.473 


546()±106 


1.2 


3.2 


-0.1 


RYVel 


28.1357 


0.424 


5124±149 


0.3 


3.5 


-0.3 


RZCMa 


4.254832 


0.084 


5800±186 


1.0 


2.7 


-0.3 


RZ Vel 


20.39824 


0.606 


4992±131 


0.5 


3.2 


-0.1 


ST Vel 


5.858425 


0.224 


5474±109 


1.1 


3.2 


-0.1 


SW Vel 


23.441 


0.583 


5150± 96 


0.5 


5.0 


-0.2 


SX Vel 


9.54993 


0.983 


6248±132 


1.3 


2.8 


-0.2 


T Vel 


4.639819 


0.135 


5800± 91 


1.8 


3.0 


-0.1 


TWCMa 


6.99507 


0.471 


5364±184 


1.0 


4.0 


-0.5 


UXCar 


3.682246 


0.577 


6002± 97 


1.7 


3.5 


-0.1 


V397 Car 


2.0635 


0.765 


5959±161 


1.5 


3.5 


-0.15 


VCar 


6.696672 


0.357 


5423±112 


1.0 


3.0 


-0.1 


V Vel 


4.371043 


0.279 


5608±121 


1.3 


3.0 


-0.3 


VXPup 


3.01087 


0.733 


5870± 97 


1.8 


3.0 


-0.15 


VY Car 


18.913762 


0.334 


6127±147 


1.5 


4.2 


-0.05 


VZ Pup 


23.171 


0.653 


5015±114 


0.2 


4.0 


-0.3 


WXPup 


8.93705 


0.497 


5263±132 


0.7 


3.0 


-0.1 



Table 3. Cepheids distances from IR photometry (d/s) and comparison with the distances computed from Berdnik ov et al. I ([2000) 
(d B ). 



Object 


d/s 


di 




kpc 


kpc 


AD Pup 


11.64 


11.03 


AH Vel 


8.59 


8.60 


APPup 


8.81 


8.75 


AQ Car 


8.20 


8.18 


AQPup 


10.15 


10.02 


AT Pup 


9.07 


8.98 


AX Vel 


8.75 


8.71 


BG Vel 


8.52 


8.51 


BN Pup 


10.64 


10.35 


DR Vel 


8.62 


8.57 


ICar 


8.40 


8.40 


MY Pup 


8.61 


8.62 


RS Pup 


9.19 


9.15 


RYCMa 


9.41 


9.26 


RYVel 


8.32 


8.31 


RZCMa 


9.99 


9.60 


RZ Vel 


8.83 


8.79 


ST Vel 


8.78 


8.68 


SW Vel 


9.06 


8.98 


SX Vel 


8.90 


8.84 


T Vel 


8.66 


8.62 


TWCMa 


10.46 


10.22 


UXCar 


8.25 


8.27 


V397 Car 


8.27 


8.27 


VCar 


8.48 


8.47 


V Vel 


8.44 


8.44 


VXPup 


9.38 




VY Car 


8.16 


8.18 


VZ Pup 


11.56 


11.35 


WXPup 


9.98 


9.80 



but the improvement was limited, since the wavelength region 
across these lines is relatively noisy in most of our spectra. 
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Table 4. Iro n, oxygen, sod i um ab unda nces, and comparison w ith p revious abundance analyses : (1) lAndrievskv et al. I ( 2002af) 
and foil. (2)lFrv & Carne vl( fT99l. (3) Luck & Lambertl dl98lh. (4) iBovarchuk & Lvubimkov.l dl984TT (5) iBarrelll dl98l. (6) 
iLuck & Lambert I (Il985l) . (7) iLuck & Lambertl (|1992|) . (8) lAndrievskv et al. I dl994l) (9) iTheveninl (1 1 998h (10) iMottini et al. I 



(2007). 



Object 


[Fe/H] 


Others 


[O/H] 


Others 


[Na/H] 


Others 




dex 


dex 


dex 


dex 


dex 


dex 


AD Pup 


-0.20±0.12 


-0.24' 


-0.09±0.04 




0.10±0.03 


0.06' 


AH Vel 


-0.06±0.05 
-0.03±O.O7 




0.14 
-0.01 




0.31±0.01 
0.25±O.O5 




APPup 


-0.18±O.O8 
-0.13±0.07 
-0.16±0.08 
-0.15±0.08 


-0.07 10 


- 
- 

-0.29 
0.08 




0.16±0.01 
0.05±0.04 
0.02±0.01 
0.03±0.06 




AQ Car 


-0.30±0.10 








0.16±0.01 




AQPup 


-0.26±0.03 


-0.14';-0.09 10 


-0.27±0.04 




0.08 


0.09 1 


AT Pup 


-0.22±0.10 




- 




0.21±0.04 




AX Vel 


-0.15±0.07 


-0.67 5 


- 




0.06±0.06 




BG Vel 


-0.10±0.08 




- 




0.16±0.06 




BNPup 


-0.03±0.10 


0.01';-0.09 10 


-0.01±0.01 


-0.06 1 


0.24±0.03 


0.09 1 


DR Vel 


-0.01±0.07 




-0.41±0.09 




0.25 




ICar 


0.11±0.07 


-0.4 4 ; -0.00 10 


-0.08 




-0.36 






0.10±0.11 


-0.3 6 


-0.37 




-0.23±0.25 






0.10±0.17 


0.02 7 


-0.12 




-0.28 






0.09±0.22 


-0.16 9 


-0.09 




0.14 




MY Pup 


-0.25±0.08 


-0.12 1 


-0.12 


-0.12 1 


0.10±0.01 


0.14 1 


RS Pup 


0.07±0.09 


Q.n u ,-0.07 3 ;-0A2 i0 


-0.02±0.22 


0.10 1 


- 


0.42 1 


RYCMa 


-0.12±0.09 
-0.20±0.09 


0.02 1 


-0.28 
- 


-0.30 1 


0.09±0.07 
-0.01±0.02 


0.16 1 


RY Vel 


-0.05±0.08 


-0.03 1 


-0.22 


-0.3 1 1 


0.31 


0.20' 


RZCMa 


-0.20±0.07 




-0.55 




-0.08±0.03 




RZ Vel 


0.05±0.10 


-0.07 1 




-0.15' 


0.29 


0.15 1 


ST Vel 


-0.14±0.10 




-0.34 




0.11 ±0.02 




SW Vel 


-0.15±0.08 


-0.07';-0.08 2 




0.18 1 




0.17 1 


SX Vel 


-0.18±0.07 


-0.03 1 


-0.31 


0.12' 


-0.01 


0.05 1 


T Vel 


-0.02±0.07 


-0.02 1 


-0.02 


0.02 1 


0.15±0.01 


0.06 1 


TW CMa 


-0.51 ±0.09 


-0.18 1 




-0.16 1 


-0.15±0.03 


-0.03 1 


UXCar 


-0.10±0.07 








0.11±0.01 




V397 Car 


-0.08±0.09 








0.24±0.06 




VCar 


-0.06±0.07 


-0.04'° 


-0.22±0.16 




0.17±0.03 




V Vel 


-0.30±0.06 








-0.34±0.17 




VXPup 


-0.15±0.12 


-0.13';-0.70 5 ; -0.39 s 


-0.10 




0.00±0.03 


0.08 1 


VY Car 


-0.06±0.06 








0.15±0.21 




VZPup 


-0.37±0.07 


-0.16 1 


-0.53 


-0.10 1 


-0.08±0.01 


o.oo 1 


WXPup 


-0.15±0.09 




-0.36 




-0.09 





4.3. Sodium 

Abundances for sodium cover the range [-0.36 dex 1-0.31 

dex]. Once again, they agree quite well with previous abun- 
dance measurements, and indeed the difference is 0.1 dex or 
less. Only two stars show a discrepancy of ~ 0.15 dex, namely 
BN Pup and RZ Vel. Most of the Cepheids in our sample are 
/Va-overabundant, a now well-known feature for Galactic super- 
giants. This TVa-overabundance is given by intermediate-mass 
MS stars, the progenitors of yellow supergiants, by the synthe- 
sis of sodium in the Ne - Na cycle in their convective cores. A 
rotationally induced turbulent and diffusive mixing brings this 
material enriched in Na into the radiative envelope and shows 
up at the stellar surface s oon after the first dredge up dSasselov I 
19861; iDenissenkov 111994 . 



4.4. Aluminium and magnesium 

Aluminium and magnesium are both under-abundant in our 
study, as expected for Galactic supergiants. As far as the Mg is 
concerned, the comparison is possible only for a few targets, and 
the abundances we obtained are quite similar to Andrievsky's 
abundances. On the other hand, the aluminum abundances 
present a larger discrepancy and current measurements are sys- 
tematically 0.1-0.2 dex lower than those by Andrievsky and 
collaborators. This systematic difference only applies to the Al 



abundances, and it could be due to a difference in the atomic 
parameters chosen for aluminium lines. 

4.5. Silicon and calcium 

These elements have, as expected, solar-like abundances. 
Current results agree quite well with the literature values, and 
indeed the discrepancies are smaller than 0.1 dex, when com- 
pared with the measurements provided by Andrievsky and col- 
laborators. In a few cases, the difference is larger, reaching ~ 0.2 
dex for Si abundances in TW CMa and for Ca abundances in 
AQ Pup, RS Pup, and VX Pup. 

5. Galactic radial abundance gradients 

iTwarog et al. I d 19971) suggests that there is a discontinuity in 
the Galactic chemical abundance gradient in the region located 
at 10-11 kpc from the Galactic centre. On the other hand, 
lAndrievskv et al. I {2004) decided to adopt a multi-modal gra- 
dient: they divided the range they studied into three regions, 
namely Zone I: 4.0-6.6 kpc, zone II: 6.6-10.6 kpc, and zone 
III: 10.6-14.6 kpc. Zone I shows a very steep gradient (-0.128 
dex kpc -1 ), the gradient becomes flatter in Zone II (-0.044 dex. 
kpc 1 ), while Zone III shows marginal evidence of a radial gra- 
dient (+0.004 dex kpc -1 ). According to these authors, there is a 
jump in the metallicity gradient in the 10-11 kpc region, where 
the iron abundance shows a rapid decrease of -0.2 dex. More 
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Tabl e 5. Aluminium, magnesium, silicon, calcium abundances and comparisons with previous abundace analyses: (1) 
lAndrievskv et al. I d2002al) . 



Object 


[Al/H] 


Others 


[Mg/H] 


Others 


[Si/H] 


Others 


[Ca/H] 


Others 




dex 


dex 


dex 


dex 


dex 


dex 


dex 


dex 


AD Pup 


-0.22±0.04 


0.00' 


-0.41 




-0.04±0.15 


-0.19' 


-0.12±0.08 


-0.17' 


AH Vel 


-0.16±0.03 
-0.13±0.04 




0.04±0.14 




0.04±0.16 
0.09±0.21 




-0.13±0.12 
-0.06±0.08 




APPup 


-0.27 
-0.17±0.06 
-0.22±0.20 
-0.24±0.03 




0.01 
-0.01 




-0.05±0.12 
0.05±0.14 
0.03±0.14 
-0.03±0.15 




-0.12±0.06 
-0.14±0.04 
-0.13±0.05 




AQ Car 


-0.19±0.05 








-0.11 ±0.24 




-0.34±0.47 




AQPup 


-0.29±0.15 


0.01 1 






-0.06±0.15 


-0.10 1 


-0.27±0.16 


-0.07 1 


AT Pup 


-0.06 








0.03 ±0.12 




-0.15±0.01 




AX Vel 


-0,35±0.04 




0.00±0.01 




-0.01±0.12 




-0.30±0.30 




BG Vel 


-0.19±0.02 




-0.34 




-0.04±0.11 




-0.16±0.09 




BNPup 


-0.10±0.03 


0.04 1 




0.21 1 


0.10±0.18 


-0.01 1 


0.05±0.01 


0.04 1 


DR Vel 


-0.08±0.07 




_ 




0.04±0.08 




-0.07 




lCar 


-0.15±0.05 
-0.43±0.20 
-0.15±0.06 
-0.13±0.06 




_ 

-0.51 
- 
- 




-0.02±0.06 
-0.16±0.15 
0.00±0.06 
0.01±0.09 




-0.16 
-0.27 
-0.08 
-0.14 




MYPup 


-0.25±0.13 


0.03 1 


-0.36±0.20 


-0.36 1 


-0.02±0.09 


-0.08 1 


-0.16±0.09 


-0.13 1 


RS Pup 


-0.08 ±0.03 


0.17 1 


- 




0.08±0.12 


0.1 1 1 


-0.02 


0.21 1 


RYCMa 


-0.10±0.14 
-0.28±0.04 




-0.05 
-0.26 




-0.04±0.11 
-0.11±0.13 


0.04 1 


-0.09±0.04 
-0.28±0.17 


-0.02 1 


RY Vel 


-0.06±0.07 


0.12 1 


_ 


-0.26' 


0.06±0.12 


0.06 1 


_ 


-0.1 1 1 


RZCMa 


-0.31 ±0.05 




-0.17±0.18 




-0.14±0.07 




-0.22±0.11 




RZ Vel 


-0.04 


0.04 1 


_ 


-0.29 1 


0.03±0.12 


0.07 1 


-0.13 


-0.07 1 


ST Vel 


-0.23±0.15 




-0.14±0.00 




-0.01±0.11 




-0.56±0.49 




SW Vel 


-0.09±0.25 


0.04 1 


- 


0.16 1 


- 


-0.03 1 


-0.27±0.21 


-0.16 1 


SX Vel 


-0.16±0.14 


0.04 1 


-0.14±0.16 


-0.20 1 


0.02±0.19 


-0.00 1 


-0.09±0.07 


-0.07 1 


T Vel 


-0.20±0.11 


0.08 1 


-0.26±0.00 


-0.24 1 


0.03±0.1 


-0.00 1 


-0.06±0.17 


-0.02 1 


TWCMa 


-0.47±0.13 




-0.45±0.00 




-0.30±0.14 


-0.10 1 


-0.24 


-0.19 1 


UXCar 


-0.30±0.01 




-0.11 ±0.1 6 




0.01±0.18 




0.00±0.01 




V397 Car 


-0.23±0.08 




-0.20 




0.00±0.28 




-0.15±0.13 




VCar 


-0.24±0.06 




-0.28±0.10 




0.01±0.25 




-0.07±0.13 




VVel 


-0.54±0.12 




-0.65±0.30 




-0.27±0.06 




-0.40±0.07 




VXPup 


-0.33±0.04 




-0.18±0.17 




-0.10±0.16 


-0.06 1 


-0.16±0.06 


-0.33' 


VYCar 


-0.16 




-0.23±0.OO 




0.00±0.24 




0.03 




VZPup 


-0.42±0.15 


-0.15' 




-0.15 1 


-0.18±0.16 


-0.09 1 


-0.35 


-0.12 1 


WXPup 


-0.29±0.04 








-0.04±0.18 




-0.21 





recent studies based o n open clusters sup port this working hy- 
pothesis. In particular, Yong et aT71 (120051) find that, instead of 
steadily decreasing with Galactocentric distance, the metallicity 
reaches a basement metallicity ([Fe/H] « -0.5 dex) beyond 10- 
12 kpc. This feature is also confirmed by the st udies of 3 field red 
giants dCarnev et al. 1120051) and of Cepheids dYong et al. f 2006) 
in the outer disk. However, in the last paper, a sub-sample of 
Cepheids reaches the basement metallicity at [Fe/H] « 0.9 
dex. By using open clusters as chemical composition tracers, 
iTwarog et al. I d2006l) found a bimodal gradient: the [Fe/H] is 
close to zero for Rc < 10 kpc, and approximately -0.3 dex for 
larger Galactocentric distances and with a marginal evidence of 
chemical composition gradients inside the two zones. 



- In the case of iron, we assigned them the value of 0.10 dex, 
since iron abundances have been determined with a suffi- 
cient number of lines, and the RMS estimates closely match 
the measurement errors. This value of 0.10 dex, greater than 
the RMS for most of our stars, is therefore slightly overesti- 
mated. 

- For the other elements, for which the number of lines used 
ranges from 1 and 5 for a large fraction of the targets, we 
assigned the value of 0.15 dex to the measurement errors. 
Such a value is very appropriate for O and Ca, and slightly 
overestimated for Na. 



The current Cepheid sample includes only a few targets in 
the crucial zone, and the use of the abundances est imated by 
lAndrievskv et al. I d2002al) and bv lMottini et al. I d2007l) does not 
allow us to provide sound constraints on the occurrence of a dis- 
continuity in the abundance gradients. Therefore, we decided to 
adopt a conservative approach by only computing linear gradi- 
ents. Abundance gradients have then been calculated with a lin- 
ear regression between abundances and Galactocentric distances 
for 5 of the 7 elements analysed in this investigation. To compute 
the goodness-of-fit, we did not adopt the errors listed in Tables 
2]and|5]because they only represent the RMS of the abundances 
when we adopt more than one line and not the real measurement 
errors that come from uncertainties in the determination of the 
atmospheric parameters. 

The measurement errors have been fixed as follow: 



To improve the sampling across the Galactic disk, we also 
included 58 C epheids from t he An drievsky's sample and 11 
Cepheids from iMottini et aTl d2007l) for which NIR photome- 
try is available in the literature. For the former sample both iron 
and heavy element abundances are available, while for the lat- 
ter only the iron abundances. The individual distances were es- 
timated using the same approach adopted in Sect. 3.7 and the 
Galactocentric distances of these additional Cepheids are listed 
in Table Q Seven of them are quite distant Cepheids located be- 
tween 10 and 1 1 .4 kpc, and a few are located toward the Galactic 
centre, thus extending the baseline we adopt to estimate the gra- 
dients. The measurement errors for these objects have been esti- 
mated by using the same assumptions we adopted for our sam- 
ple. 
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5. 1 . Iron Galactic radial gradient 

We estimated at first the linear Galactic radial gradients for 
the current Cepheid sample. The slopes, their 1 — <x errors, 
and the correlation coefficients are listed in columns 2,3, and 
4 of Table [6] Data plotted in the top panel of Fig. [7] show 
that the slope is -0.061 dex kpc" 1 in the region ranging from 
ss 8 to 12 kpc. This seems to be a relatively high value 
when comp ared with the slop e s deter mined in a similar radial 
interval by lAndrievskv et al. I d2002 aj). who found -0.029 dex 
kpc" 1 under the hypothe sis of a single-zone gradient, and by 
AndrievskxitaD (12004 : -0.044 dex. kpc 1 for a multimodal 
gradient. We note in passing that the correlation coefficient is 
relatively small (r =0.44), but quite s i milar to the correlation co- 
efficient found bv lAndrievskv et al. I d2002al) (r=0.47). 

In order to ascertain on a quantitative basis whether the cur- 
rent slope evaluation is affected by the size of the sample, we also 
included the Cepheids from Andriev sky's sample (29 Cepheids) 
and from the iMottini et aTl d2007l) sample (5 Cepheids) with 
Galactocentric distances ranging from ~ 8 to ~ 12 Kpc. This 
means that we doubled the sample size. The results of the lin- 
ear regression are listed in columns 5, 6, and 7 of Table [6] 
Interestingly enough, data plotted in the top panel of Fig.[8]show 
that the new slope, within the I — cr error, agrees quite well (- 
0.056 vs -0.061 dex kpc 1 ) with the slope only based on the cur- 
rent Cepheid sample. The same outcome applies to the correla- 
tion coefficient. 

Futhermore, we also decided to check whether current slope 
evaluations are affected by the Galactocentric radial distances 
covered by current Cepheid sample. Therefore, we included 
the entire Andrievs ky sample (58 Cepheids) and the entire 
Mottini et al. (2007) sample (11 Cepheids). This means that we 
increased the size of the sample by more than a factor of three 
and by more than 40% the range of Galactocentric distances. The 
results of the linear regression over the entire sample are listed 
in Table [8] The top panel of Fig. [9] once again shows that the 
slope based on the entire data sample agree, within the errors, 
with those based on a smaller Cepheid sample. However, the 
1 - cr error is now a factor of two lower (0.008 vs 0.19); and the 
correlation coefficient increased by almost the 50%. These find- 
ings bring out several circumstantial evidence: i) - the slope of 
the iron radial gradient appears marginally affected by the uncer- 
tainties in radial distances, and indeed the current spread could 
be due to uncertainties in iron abundances; 
ii) - The correlation coefficient depends on the intrinsic accuracy 
of individual distances, on the range of Galactocentric distances 
covered by the sample, and to a lesser extent on the spread in 
metal abundances. 

Hi) - The slope based over the entire Cepheid sample is 
marginally affected by the few objects located toward the 
Galactric centre (d < 6 Kpc) and toward the outer reaches of 
the Galactic disc (d > 1 1 Kpc. In fact, we performed a linear 
regression by removing these objects and both the slope and the 
1 - cr error are marginally affected. 

iv) - The current iron radial gradient agrees quite well wi th pre- 
vious determinations. In particular, Twaro g et al. I (1 19971) com- 
puted a single-zone gradient using open clusters covering a 
broad radial interval from 6 to 16 kpc, and they found -0.067 
dex kpc 1 . The gradient changes to -0.071 dex kpc 1 if the clus- 
ter BE 21 is included in the sample with a high metallicity value 
and to -0.076 dex kpc -1 if it is included with a low metallic- 
ity value^Awn^similar value (-0.06 dex kpc -1 ) was also found 
bv lFriel et al.1 d2002l) using the same chemical composition trac- 
ers. By adopting classical Cepheids, Caput oet al. I d200ll) found 



a slope of -0.05 dex kpc 1 , which is slightly smaller than the 
current estimate, but still within the error bars. Our result is also 
fully consistent wit h one of the las t gradi ent values determined 
with Cepheids by iKovtyukh et al. I d2005l) . who added data for 
16 distant Cepheids to the Andrievsky's sample. By choosing 
a single-zone gradient in the 4-15.5 kpc range, they obtained a 
gradient of -0.06 dex kpc 1 . Finally, we can compare ou r slope 
with the very recent determination bv lLuck et al. I ((2006). They 
found an iron gradient of -0.068 dex kpc 1 , a value that is in re- 
markable agreement with current estimate. The radial gradients 
based either on classical Cepheids o r on open clusters disagree 
only w ith the estimates provided by iKilian-Montenbruck et al. I 
dl994h . who found a gradient of -0.003+0.020 dex kpc" 1 , using 
B-type stars in the 6-15kpc region. 



5.2. Galactic radial gradient for other elements 

Galactic radial abundances gradients for the other elements 
in the 8-12 kpc region, computed either from current sam- 
ple or from the sample implemented with Cepheids from the 
Andrievsky's sample are once again in good agreement. The dis- 
crepancies between the two different slopes are on average of the 
order of 0.01 dex kpc" 1 . The only exception is for Na and O, for 
which the discrepancy is ~ 0.015 dex kpc" 1 , but still within the 
error bars. 

The radial gradients for the other elements in the radial 
range from 8 to 12 kpc present larger uncertainties when com- 
pared with the iron one . A similar trend was also found by 
lAndrievskv et al. I d2002al) . Once again, curren t slopes are higher 
when c ompared with those estimated by lAndrievskv et al. I 
(l2002al) . In particular, for oxygen they found -0.022+0.009 
dex kpc" 1 , while we find a steeper slope -0.041+0.034 dex 
kpc" 1 . The same outcome applies to Na for which they found 
-0.023+0.006 dex kpc" 1 , while we find -0.042+0.029 dex kpc" 1 . 
On the other hand, the slopes for Si and Ca agree very 
well, and indeed the slopes are -0.030+0.004 dex kpc" 1 and - 
0.021+0.006 dex kpc" 1 , while we find -0.031+0.029 dex kpc" 1 
and -0.014+0.029 dex kpc" 1 . 

If we compute the radial gradients using the entire Cepheid 
sample, they appear to be similar to the iron gradient and range 
from -0.06 to -0.07 dex kpc" 1 . The gradient appears to be flat- 
ter across the solar circle. This behavio ur was expect e d, be cause 
it has already been detec t ed, firs t by [Twarog et al. I d 1997b but 
also bv lAndrievskv et al7l d2002al) . This circumstantial evidence 
supported the use of a multi-zonal gradient. However, before 
firm conclusions can be derived, a large sample of homogenous 
chemical abundances for iron and heavy elements is required. 
As far as oxygen is concerned, the co mparison with previ- 



ous st udi es shows a good agr ee ment with Smartt & Rolleston 
(119971) . iGummersbach et al. I d!998l) . and iRolleston etaL 
(I2OO0I) . who found slopes of -0.07+0.01 dex kpc" 1 , - 
0.067+0.024 dex kpc" 1 , and -0.067+0.008 dex kpc" 1 , using 
B-type stars as c h emical composition tracers. Futhermore, 
iMaciel&Ouirezal dl998l) f ound -0.058+0.007 de x kpc" 1 
using planetary nebulae and Afflerba ch et al. I (1 19971) found - 
0.064+0.009 dex kpc" 1 from HII regions. However, current 
results do not sup port the flat slope s for t he oxygen radial 
gradient fou nd by [F itzsimmons et al. I (fl992l) (-0.03+0.02 dex 
kpc" 1 ) and iKilian-Montenbruck et all d 19941) (-0.021 ±0.0 12 
dex kpc" 1 ) using B-tvpe stars or by iDeharveng et al. I ((2000) 
(-0.0391+0.005 dex kpc" 1 ) using HII regions. Obviously 
they also do not ag ree with the null radial gradient found by 
iKaufer et al. I (1 19941) (-0.000+0.009 dex kpc 1 ) using B-type 
stars. 



Lemasle et al.: Detailed chemical composition of Galactic Cepheids. 9 

Table 6. Galactic radial abundance gradient in the 8-12 kpc region. The left part of the tables are the valu es for our sample, the right 
part includes data from Andr ievskv et al. I d2002al) and subsequent papers and from lMottini et al.1 (120071) . 



Element 


slope 
dex kpc -1 


l-o- 
dex kpc -1 


Correlation coeff. 


slope 
dex kpc -1 


1-cr 
dex kpc -1 


CoiTelation coeff. 


Fe 


-0.061 


0.019 


0.441 


-0.056 


0.012 


0.410 





-0.041 


0.034 


0.226 


-0.051 


0.022 


0.305 


Na 


-0.042 


0.029 


0.269 


-0.027 


0.020 


0.199 


Si 


-0.031 


0.029 


0.319 


-0.044 


0.020 


0.433 


Ca 


-0.014 


0.029 


0.101 


-0.024 


0.020 


0.160 



Tabl e 7. Galactoce ntric distances for Cepheids in And rievskv et al. I (l2002al) and following (3 first couples of columns) and in 
Mott ini et al. I d2007l) (last couple of columns) for which distances based on IR photometry are available. 



Object 


Distance 


Object 


Distance 


Object 


Distance 


Object 


Distance 




kpc 




kpc 




kpc 




kpc 


SZ Aql 


6.96 


XLac 


9.01 


AV Sgr 


6.37 


U Car 


8.11 


TT Aql 


7.7 


Y Lac 


9.08 


BB Sgr 


7.7 


WZCar 


8.14 


FM Aql 


7.86 


Z Lac 


9.2 


V350 Sgr 


7.61 


VWCen 


6.83 


FN Aql 


7.41 


BG Lac 


8.77 


RV Sco 


7.74 


XX Cen 


7.52 


V496 Aql 


7.63 


TMon 


9.75 


RYSco 


7.16 


KNCen 


6.92 


V1162Aql 


7.43 


SV Mon 


10.83 


KQ Sco 


6.26 


GH Lup 


7.56 


Eta Aql 


8.3 


CV Mon 


10.09 


V500 Sco 


7.08 


SMus 


8.12 


V340 Ara 


4.99 


S Nor 


7.73 


SS Set 


7.61 


UU Mus 


7.6 


RT Aur 


8.94 


V340 Nor 


6.9 


UZSct 


5.69 


U NOI- 


7.34 


SU Cas 


8.71 


YOph 


7.93 


EWSct 


8.14 


LS Pup 


11.39 


V Cen 


7.99 


BFOph 


7.68 


V367 Set 


6.81 


EU Tau 


8.87 


Del Cep 


8.57 


GQ Ori 


10.84 


BQ Ser 


7.7 






BG Cm 


8.3 


XPup 


10.5 


STTau 


9.53 






XCyg 


8.32 


SSge 


8.13 


SZ Tau 


8.96 






VZCyg 


8.75 


USgr 


7.89 


S Vul 


7.6 






Beta Dor 


8.5 


WSgr 


8.08 


T Vul 


8.34 






WGem 


9.45 


YSgr 


8.02 


U Vul 


8.16 






RZ Gem 


10.66 


VY Sgr 


6.29 


SV Vul 


7.8 






AD Gem 


11.39 


WZ Sgr 


6.71 










Zeta Gem 


8.86 


APSgr 


7.68 











Tabl e 8. Galactic radia l abundance gradient in the 5-12 kpc region . Our sample has been completed with Cepheids from 
lAndrievskv et al. I d2002al) and following and from lMottini et al.1 ((2007) for which a homogeneous determination of NIR distances 
was possible. 



Element 


slope 


l-o- 


Correlation coeff. 




dex kpc -1 


dex kpc -1 




Fe 


-0.070 


0.008 


0.633 





-0.065 


0.013 


0.471 


Na 


-0.071 


0.013 


0.524 


Si 


-0.063 


0.012 


0.666 


Ca 


-0.062 


0.012 


0.484 



For silicon , our r esults are in good agreement with 
iRolleston et ail d2000l) . who found a slope of -0.06+0.01 
dex kpc" 1 , but agree n e ither with the flatter slopes by 
iKilian-Montenbruck et al. I dl994l) (+0.000+0.018 dex kpc ' 1 ) 
nor with the steeper slopes by fcummers bach etaTI[T998l) (- 
0.107+0.028 dex kpc" 1 ). Both these authors were using B-type 
stars. 

However, in order to ascertain whether the difference in the 
different slope estimates are dominated by deceptive empirical 
uncertainties or they might be intrinsic, a detailed investigation 
of stellar ages and radial distribution of the different chemical 
tracers is required. 



6. Conclusions 

We have determined the abundances of iron and six light met- 
als (O, Na, Mg, Al, Si, Ca) in 30 Galactic Cepheids. These 
abundances were used to compute a slope for the Galactic ra- 
dial abundance gradient. For the iron Galactic gradient in the 
radial distance range from 8 to 12 Kpc we found slopes of 
-0.061 dex kpc" 1 and of -0.056 dex kpc" 1 if we include 69 
Cepheids with homogeneous distance determinations. The cur- 



re nt slope is steeper when c ompared with the slope estimated 
by [Aridrievskv et al. I (12002a) for the central region of the multi- 
zonal gradient. On the other hand, the current gradient computed 
over a more extended baseline (-0.070+0.008 dex kpc" 1 ) is fully 
consistent with other studies based either on open clusters or on 
Cepheids. Current findings do not allow us to provide firm con- 
straints on the multi-modal model, due to the limited number of 
targets in the crucial zone. The hypothesis of a multi-modal gra- 
dient deserves further investigations, mainly to fill the transition 
region across 10-11 kpc. 
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Fig. 1. Observed curve of growth for Fe I in APPup. The full line is the theoretical curve of growth for a typical line (/i=5000A, 
Xex=$)- The atmospheric parameters adopted for this star are T c //=6180 K, log g=2.2, v ( =4.0 km/s, and [Fe/H]=-0.2 dex. The 
dashed line on the left shows the [Fe/H]=0 location. 




lo< 3 («© * * r ) 

Fig. 2. Observed curve of growth for Fe II in APPup. The full line is the theoretical curve of growth for a typical line (/l=5000A, 
Xex=3)- The atmospheric parameters adopted for this star are T e //=6180 K, log g-2.2, v,=4.0 km/s, and [Fe/H]=-0.2 dex. The 
dashed line on the left shows the [Fe/H]=0 location. 
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Fig. 3. Observed curve of growth for Fe I in TW CMa. The full line represents the theoretical curve of growth for a typical line 
(^=5000A 

■>Xex— 3)- The atmospheric parameters adopted for this star are T e ff— 5364 K, log g-1.0, v r — 4.0 km/s, and [Fe/H]— 0.5 
dex. The dashed line on the left shows the [Fe/H]=0 location. 
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Fig. 4. Observed curve of growth for Fe II in TW CMa. The full line represents the theoretical curve of growth for a typical line 
(^=5000A 

■>Xex— 3)- The atmospheric parameters adopted for this star are T e ff— 5364 K, log g-1.0, v r — 4.0 km/s, and [Fe/H]— 0.5 
dex. The dashed line on the left shows the [Fe/H]=0 location. 
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Fig. 5. Observed curve of growth for Fe I in VZPup. The full line represents the theoretical curve of growth for a typical line 
(i=5000A 

■>Xex— 3); The atmospheric parameters adopted for this star are T e yy— 5015 K, log g— 0.2, v ( -4.0 km/s, and [Fe/H]— 0.3 
dex. The dashed line on the left shows the [Fe/H]=0 location. 
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Fig. 6. Observed curve of growth for Fe II in VZPup. The full line represents the theoretical curve of growth for a typical line 
(^=5000A 

■>Xex— 3)- The atmospheric parameters adopted for this star are T e ff— 5015 K, log g-0.2, v,-4.0 km/s, and [Fe/H]— 0.3 
dex. The dashed line on the left shows the [Fe/H]=0 location. 
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Fig. 7. Galactic radial abundance gradients in the [8-12] kpc region based on our Cepheid sample. The slopes, the 1 — cr errors, and 
the correlation coefficients are also labelled. The solid line shows the linear regression. 
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Fig. 8. Galactic radial abundance gradients in the T8-121 kpc region b ased on our Cepheid sample (crosses) and on data from 
Andrievsky and collaborators (triangles), and from Mottini et al. I d2007l) (squares). 
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Fig. 9. Same as Fig. 8, but for Cepheids with Galactocentric distances ranging from 5 to 12 kpc. 



